Periodic Anderson model with electron-phonon correlated conduction band 
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This paper reports dynamical mean field calculations for the periodic Anderson model in which 
the conduction band is coupled to phonons. Motivated in part by recent attention to the role 
of phonons in the 7-a transition in Ce, this model yields a rich and unexpected phase diagram 
which is of intrinsic interest. Specifically, above a critical value of the electron-phonon interaction, 
a first order transition with two coexisting phases develops in the temperature-hybridization plane, 
which terminates at a second order critical point. The coexisting phases display the familiar Kondo 
screened and local moment character, yet they also exhibit pronounced polaronic and bipolaronic 
properties, respectively. 

PACS numbers: 71.27. +a, 71.10.Fd, 71.38.-k 



The Periodic Anderson Model (PAM) and its impurity 
variant have played pivotal roles in elucidating the na- 
ture of Kondo screening as the techniques of many-body 
theory have improved P, Q • Perhaps its most noted ap- 
plication has been the Kondo Volume Collapse scenario 
for understanding the unique isostructural 7-a transition 
in Ce, with its very large 15% volume change The 
relative merits of this perspective versus the Mott tran- 
sition scenario [j| are still under debate, although both 
focus on critical 4/-electron correlation effects, and the 
finite temperature predictions are rather similar Q . The 
PAM exhibits a smooth crossover from a local moment 
region with Curie Weiss susceptibility (7-like) to a re- 
gion with Kondo screened 4/ moments and a paramag- 
netic susceptibility (a- like), as a function of increasing 
hybridization between the 4/ and valence electrons 0- 
Although the PAM also predicts a first order transition 
given proper consideration of the Maxwell construction of 
the free energy versus volume curves [3| , it requires modi- 
fications like some /— / hybridization |8j so as to display a 
first order phase transition with two coexisting phases at 
the same hybridization in the temperature-hybridization 
plane. 

Over the past decade attention in the Ce literature 
has shifted to an appreciation that a significant fraction 
of the total entrop y ch ange across the transition may be 
due to phonons [9l4l5l|. However, studies focusin g on the 
effect of phonons on the PAM are very limited p3S|. 
Prior studies either are constrained to ground state cal- 
culation or do not explore possible phase transitions in 
detail. To this end, we are motivated here to consider 
the PAM with Holstein phonons 21-24]]. Since the cou- 



pling of phonons to the /-electrons can lead to loss of 
local moments via electron condensation, we have cho- 
sen to couple the phonons to the conduction electrons 
in the present work. We find that the electron-phonon 
interaction above a critical strength induces a first order 
transition in the temperature-hybridization plane for the 
PAM-Holstein model. Strikingly the electron-phonon in- 



teraction also creates polaronic behavior in the Kondo 
screened phase and bipolaronic behavior in the local mo- 
ment phase. This intriguing phase diagram is explored 
in the remainder of the present paper. 

The Hamiltonian of the PAM-Holstein model is: 

H = H + Hjj + £f e _ p h, (1) 
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where c^ a , c\ a (fi,o-, f} a ) arc the creation and annihi- 



lation operators of the conduction (/-level) at site i and 
spin a; n\ a = \ a c i>a and n{ CT = fj^fi,a represent the 
occupation of the c and /-electrons, respectively; t is the 
nearest neighbor hopping; e/, the on-site energy of the 
/-level; V, the hybridization between conduction and lo- 
calized electrons; the on-site Hubbard interaction is U; g 
is the electron-phonon coupling; JQ, the lattice displace- 
ment at site i, and Pi its conjugate momentum. 

We use the Dynamical Mean Field Theory (DMFT) 
25j on a hypercubic lattice in infinite dimensions with 

Gaussian density of states Die) — —= — e~( w ) . We 
J v ; sf^W 

set the bandwidth W = 1.0 as the unit of energy. In 
Ce the Fermi energy is about 6000 K and the Debye fre- 
quency is 110-160K [j| HH, therefore we set the phonon 
frequency ujq = 0.01 at 1% of bandwidth. The total 
electronic density is fixed at n = 1.8 by tuning the chem- 
ical potential at each iteration of the DMFT cycle. The 
Hubbard interaction is U = 4.0, and we adjust e/ so 
that ri/ = latT = 0.1to ensure that a local moment 
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is present at large temperatures. Therefore the data we 
show are for nf ~ 1.0 and n c ~ 0.8. The Continuous 
Time Quantum Monte Carlo (CTQMC) [H], generalized 
for electron-phonon coupling [27J , is employed as the im- 
purity solver. 

Fig. [T] displays the local hybridization factor T = 
( c o fo + h.c.) (here denotes the impurity site) as a func- 
tion of V for g 2 /2k = 1.0 and different values of inverse 
temperature, /?. As the temperature decreases, the slope 
of the r vs. V curve becomes progressively larger, which 
indicates the system is approaching a critical point. In- 
terestingly, the curves approximately cross for a critical 
hybridization of V c ~ 0.96. The inset of Fig. [1] shows T 
vs. V at g 2 /2k = 0.49. Notice that for this value of the 
coupling the slope docs not become steeper as the tem- 
perature decreases, and the crossing disappears. This 
indicates that the corresponding susceptibility reaches a 
plateau as a function of the temperature and the criti- 
cal behavior is lost. We believe that g 2 /2k = 0.49 is a 
lower bound for the critical value of the electron-phonon 
coupling. 




FIG. 1: (color online) Isothermal scan of the hybridization 
factor r = (cq fo + h.c.\ as a function of V for g 2 /2k = 1.0. 
r increases monotonically with V. As the temperature de- 
creases, F vs. V becomes steeper with a diverging slope near 
V c ~ 0.96. Inset: The isothermal scan of the hybridization 
factor F as a function of V at g 2 /2k — 0.49. Notice that the 
critical behavior has disappeared. 

When the temperature is further decreased to T = 
0.0167 ((3 = 60), r vs. V displays a hysteresis loop as 
shown in Fig. [2J The red line is obtained by starting at 
the large V side (V = 1.2), and using the output self- 
energy to initiate the simulation for the next smaller V. 
On the other hand, we obtain the black line by starting 
at V = 0.8 and using the output self-energy as the input 
for the next larger value of V. The coexistence of two 
solutions for the same value of V at T = 0.0167 is a direct 
evidence of a first order phase transition. The absence of 
such a hysteresis at higher temperatures indicates that 
the first order transition ends at a second order terminus 
(Vc, T c ). 

For the same parameters, V = 0.96, g 2 /2k = 1.0, u; = 
0.01, and U = 4.0, we also perform a series of isothermal 



scans on the chemical potential to study the relationship 
between the total electron density n = n c + Uf and the 
chemical potential [i. As long as the temperature is not 

din 

below T = 0.0167, the compressibility — — shows no ten- 

dfi 

dency to diverge. This indicates the phase transition here 
is not compressibility driven. 




FIG. 2: (color online) Hysteresis of T vs. V for T = 0.0167, 
g 2 /2k = 1.0. The black line represents the small V branch of 
the hysteresis for which the self-energy of the previous sim- 
ulation is used to initiate the calculation for the next larger 
value of V . While the red line represents the large V branch 
where starting with V = 1.2 we use the output of the previous 
simulation to initiate the computation at the next lower value 
of V. Inset: T as a function of the DMFT iteration number n 
for V = 0.96, T = 0.0167. The black (red) symbols represent 
the small (large) V branches. 

In Fig. [3] we show the temperature times the local f- 
orbital spin susceptibility, T ■ x{^, versus temperature. 
As T approaches zero T ■ xl^ is roughly constant for 

V = 0.8, while it goes to zero for V = 1.2. This indicates 
that at V = 0.8 the /-electrons display a robust local 
moment and paramagnetic local susceptibility with 1/T 
dependence, while at V — 1.2 the / local moments are 
quenched. The inset of Fig. [3] shows the /-orbital density 
of states (DOS) at T = 0.01. Notice that at V = 0.8 
there is a gap across the Fermi level, while at V = 1.2 
a Kondo resonance peak appears. The screening of the 
local moment in the large V region is a consequence of 
the singlet formation between c and /-electrons. 

The main panel of Fig. [5] shows the occupancy dis- 
tribution histogram of the c-electrons, P(n c ), at T = 
0.0167. P(n c ) has been used to illustrate bipolaron for- 
mation [27J. At V = 0.8 the c-orbital electrons are in 
a bipolaronic state, which is characterized by the oscil- 
lation between zero and double occupancy. While for 

V = 1.2, the c-electrons are in a polaronic state, where 
the occupancy oscillates between zero and one. For the 
PAM, without electron-phonon coupling, the structure of 
P(n c ) is totally different. Here there is only one peak at 
roughly the c-electron filling n c = 0.8, and P(n c ) quickly 
decays to zero for n c away from this filling. In the in- 
set, the quasi-particle fraction Z is plotted as a function 
of temperature. The quasi-particle fraction is calculated 
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FIG. 3: (color online) Temperature times the local /-orbital 
spin susceptibility, T ■ xV\ as a function of temperature for 
g 2 /2k = 1.0. For V = 0.8 (black line), T ■ xl f approaches a 
constant value as T — > indicating an unscreened moment. 
For V = 1.2 (red line), T ■ xl^ converges to zero indicating 
the local moment is screened. Inset: The /-electron DOS at 
T = 0.01. The Kondo peak found for V = 1.2 (red line), but 
absent for V = 0.8 (black line) is consistent with the screening 
and unscreening scenarios in the main panel. 

for the lower quasiparticle band at the Fermi level us- 
ing a generalization of the single band formulation pjf . 
The main component of this approach is to make the re- 
di?eE(w) ImE(wrT) 



placement 



which becomes 
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exact at zero temperature. As T — > 0, Z goes to zero for 
V = 0.8 indicating non-Fermi liquid behavior, while it 
converges to a finite value for V = 1.2, the signature of 
Fermi liquid formation. 
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FIG. 4: (color online) Occupancy distribution histogram of c- 
orbital P(n c ) for V = 0.8 (black line) and 1.2 (red line), T = 
0.0167 and g 2 /2k = 1.0. For comparison, P(n c ) of the PAM 
without electron-phonon coupling is plotted as well: V = 0.8 
(blue line) and V = 1.2 (green line). Inset: the quasiparticle 
Z factor as a function of temperature for V = 0.8 (black line) 
and 1.2 (red line). 

We find that this Kondo singlet to local moment phase 
transition remains for a large range of parameters, like 
adjusting the total filling to n = 1.6, changing the Hub- 
bard interaction to U = 3.8 and increasing the phonon 
frequency to ujq = 0.02 and u>o = 0.05, while keeping 
g 2 /2k fixed. For these different parameters, we find that 



the isothermal T vs. V curves still cross and their slopes 
diverge at a critical value of the hybridization, V c , as the 
temperature is decreased. We also find that V c changes 
roughly linearly with g 2 /2k. 

In Fig. Eta) the time integrated local /-orbital spin- 
spin correlation function, is plotted as a function 
of temperature for V = 1.1,1.2 and 1.3. We identify 
the Kondo scale Tk as the energy where xl ^ a ^ s to 
around half of its low temperature value. We find that 
Tk changes very little as V increases, so the line V vs. 
Tk should have a large slope. Fig. [S^b) shows the time 
integrated local c-orbital spin-spin correlation function. 
Xs C vs. V, at different temperatures, where large values 
reflect the c-electron spin degeneracy in the polaronic 
state in contrast to the small susceptibility for the spin- 
less bipolarons. For V < 0.96 the curves almost overlap 
for all T < 0.1. In fact, the corresponding c-electron oc- 
cupancy histograms (not shown) show an obvious bipola- 
ronic double peak feature even at relatively high temper- 
atures like T — 0.1. If we define T* as the energy where 
bipolaron formation begins, then the line T* vs. V must 
be nearly horizontal. 
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FIG. 5: (color online) Panel (a) shows the /-orbital time 
integrated local spin-spin correlation, xl \ as a function of 
temperature for V = 1.1, 1.2 and 1.3. Panel (b) shows the 
c-orbital time integrated local spin-spin correlation function, 
XT, as a function of V for different temperatures. 

We have also calculated the renormalized phonon fre- 
quency. At T = 0.025 it is roughly constant for hy- 
bridization V > 0.96; however, it drops precipitously for 
V < 0.96, decreasing by half when V = 0.8. This be- 
havior softens with increasing temperature, e.g., a more 
gradual decrease begins for V < 1.2 at T = 0.1. This 
indicates an important temperature dependence of the 
phonons properties. Indeed the analysis in [l3j for Ce 
found that the temperature dependence of the phonons 
was a critical factor for obtaining a significant phonon 
contribution to the entropy change across the 7-cv tran- 
sition [l-Ell]. 

Fig. |5] is a schematic summary of our findings. Two 
phases, local moment-bipolaron and Kondo singlet- 
polaron, are separated by a first order transition line, 
which terminates at a second order critical point {V c , T c ). 
The positive slope of the V vs. T first order transition 
line is a consequence of a Clausius-Clapeyron-like relation 
where hybridization V is the intensive analog of pressure. 
There is no broken symmetry between these two phases 



4 



V 



Kondo Singlet 
Polaron 



Local Moment 
Bipolaron 



Tk 



Critical Point 



FIG. 6: (color online) Schematic V vs. T phase diagram. 
The solid black line represents the first order phase transition 
which separate the local moment-bipolaron phase for small V 
from the Kondo singlet-polaron phase for large V. This first 
order phase terminates at a second order critical point. The 
red dashed line coming out of the critical point represents 
the Kondo scale Tk and the green dashed line the bipolaron 
energy scale T* . 

as we can move adiabatically from one to another by 
wandering around the critical point. Both phases are de- 
stroyed by increasing the temperature. In order to have 
such a first order phase transition, the electron-phonon 
coupling on the c-band must be larger than a certain 
critical value. The fact that the critical temperature is 
a function of electron-phonon coupling implies that the 
critical point touches zero temperature at some g c , where 
the first order phase transition becomes a quantum phase 
transition tuned by V. 

In conclusion, when the conduction band of the peri- 
odic Anderson model is coupled to phonons, one obtains 
a rich and unexpected phase diagram. Above a criti- 
cal strength of the electron-phonon coupling a first or- 
der transition with two coexisting phases develops in the 
temperature-hybridization plane. This transition termi- 
nates at a second order critical point. These coexisting 
phases correspond to the familiar Kondo screened and 
local moment regions of the PAM, yet, they additionally 
exhibit pronounced polaronic and bipolaronic behavior, 
respectively. While the PAM and its impurity variant 
have been paradigms for the 7-0 transition in Ce, ad- 
ditional electronic bands not considered here might be 
needed in a generalization of the present PAM-Holstein 
model to more completely explain the volume collapse. 
Nonetheless, the present simple model illustrates a fun- 
damental principle relevant to the electron-phonon inter- 
action in Ce. The Kondo temperature, which measures 
the critical energy scale of hybridization between 4/ and 
valence electrons, has a roughly exponential volume de- 
pendence leading to an order of magnitude increase from 
7- to a-Ce [1, U This scale is comparable to that of the 



lattice vibrations (Debye temperature) [15| only in the 7 
phase of Ce, and so it is no accident that the present work 
finds the most dramatic manifestations of the electron- 
phonon interaction, a bipolaronic state with significant 
phonon softening, in this local moment region. 
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